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Abstract 



B 

^ . We study the coupled particle and energy transport in a prototype model of 

G I interacting one- dimensional system: the disordered hard-point gas, for which 

Q I numerical data suggest that the thermoelectric figure of merit ZT diverges 

with the system size. This result is explained in terms of a microscopic 
mechanism, namely the local equilibrium is characterized by the emergence 
of a broad stationary "modified Maxwell-Boltzmann velocity distribution", 

T^ ! of width much larger than the mean velocity of the particle flow. 
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1. Introduction 

^ . Thermoelectricity is an old field: The Seebeck effect, that is, the con- 

H i version of temperature differences into electricity, was discovered in 1821. 

However, a strong interest in termoelectric phenomena arose only in the 
1950's, when loffe discovered that doped semiconductors exhibited much 
larger thermoelectric effect than did other materials. He also proposed that 
semiconductors could be used to build solid-state home refrigerators. Such re- 
frigerators would be long-lived, silent, maintenace-free, and environmentally 
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benign. loffe's suggestion initiated an intense research activity in semicon- 
ductors physics pj, y, |3|. However, in spite of all efforts and consideration of 
all type of semiconductors, thermoelectric refrigerators have still poor effi- 
ciencies compared to compressor-based refrigerators. Today, thermoelectric 
devices are mainly used in situations in which reliability and quiet opera- 
tion are more important than the cost. Applications include equipments in 
medical applications, space probes, etc. 

In the last decade there has been an increasing pressure to find better 
thermoelectric materials with higher efficiency. The reason is the strong 
environmental concern about chlorofluorocarbons used in most compressor- 
based refrigerators. Also the possibility to generate electric power from waste 
heat using thermoelectric effect is becoming more and more interesting ll|, |2|, 

The thermodynamic efficiency can be conveniently written in terms of 
the so-called figure of merit ZT = {aS"^ / k)T, where a is the electric con- 
ductivity, S the thermoelectric power (Seebeck coefficient), k the thermal 
conductivity, and T the temperature. Ideal Carnot efficiency is recovered in 
the limit ZT — )■ oo. In spite of the worldwide research efforts for identifying 
thermoelectric materials with high ZT values, so far the best thermoelec- 
tric materials are characterized by values of ZT ~ 1, at room temperature. 
Values ZT > 3 are considered to be essential for thermoelectric devices to 
compete in efficiency with mechanical power generation and refrigeration. 

The challenge lies in engineering a material for which the values of a, S, 
and K, can be controlled in order to optimize thermoelectric efficiency. The 
problem is that the different transport coefficients are interdependent, thus 
making optimization extremely difficult. On the other hand, thermodynam- 
ics does not impose any upper bound on ZT, so that efficient thermoelectric 
devices could in principle be engineered. The present understanding of the 
possible microscopic mechanisms leading to an increase of ZT is quite lim- 
ited, with few exceptions. Notably, Refs. [sl, |6[ showed that the optimal 
density of states in a thermoelectric material is a delta function. Such sharp 
energy filtering allows to reach, in principle, the Carnot efficiency. 

Here we consider the problem of increasing thermoelectric efficiency from 
a new perspective, that is, we pursue a dynamical system approach. Under- 
standing from first principles and from nonlinear dynamics simulations the 
microscopic mechanisms that can be implemented to control the heat flow |7| 
might prove useful not only for thermoelectric phenomena but also for the 
design and engineering of thermal diodes and transistors. In this paper, our 



plan is to compute transport coefficients and thermoelectric efficiency from 
ffist principles, namely from the underlying microscopic dynamical processes 
which are known to be predominantly nonlinear in nature. In a previous 
work J8|, the thermoelectric problem has been investigated by numerical 
solution of the microscopic equations of motion. Inspired by the kinetic 
theory of ergodic gases and chaotic billiards, a simple microscopic mecha- 
nism for increasing thermoelectric efficiency was proposed. More precisely, 
the cross transport of particles and energy in open classical ergodic billiards 
was considered. It has been shown that, in the linear response regime, the 
thermoelectric efficiency can approach Carnot efficiency for sufficiently com- 
plex charge carrier molecules. Indeed, the figure of merit has been found 
to be a growing function of the number di^t of internal degrees of freedom, 
ZT = {d + 1 + (iint)/2, where d is the geometric dimension. 

In spite of the abstract nature of the model, the above paper opens the 
possibility for a theoretical understanding of the basic microscopic require- 
ments that a classical dynamical system must fulffil in order to lead to a high 
thermoelectric figure of merit. In particular, the question arises whether 
inter-particle interaction might increase the effective number of degrees of 
freedom, thus leading to a higher figure of merit than in the noninteracting 
idealized d dimensional gas, where ZT = {d + l)/2. Alnog these lines a de- 
tailed numerical study of the cross heat and particle transport has been per- 
formed for an open one-dimensional disordered hard-point gas J9|. It has been 
found that ZT diverges as a power-law in thermodynamic limit, ZT oc A^*, 
where N is the average number of particles in the system and b ^ 0.79. Even 
though the above result could be, in principle, very interesting, no indication 
was given concerning the microscopic mechanism which is responsible for the 
increase of ZT. On the other hand a theoretical understanding is needed in 
order to obtain useful hints for increasing thermodynamic efficiency in more 
realistic models. 

In this paper, we propose a mechanism which explains the large thermo- 
electric quality factor ZT numerically observed in Ref. |9|. This mechanism 
requires local equilibrium, as naturally expected in systems with the mix- 
ing property, and the emergence, in the linear response regime, of an out- 
of-equilibrium "modified Maxwell-Boltzmann velocity distribution" of width 
much larger than the mean velocity of the particle flow. Such broad distribu- 
tion limit is opposite to the limit of peaked distribution, corresponding to the 
delta-like energy filtering put forward in Refs. p, |6|. We provide numerical 
evidence supporting the effectiveness of the broad-distribution mechanism in 



the hard-point gas modeL 

Our paper is organized as follows. Sees. HI El and H] review introductory 
material on coupled particle and energy transport, modeling stochastic baths, 
and thermoelectric efficiency of the one-dimensional ideal non- interacting gas. 
In Sec. O we present numerical results for thermodynamic transport coeffi- 
cients for the disordered hard-point gas model. Finally, the obtained numer- 
ical results are explained in terms of a mechanism based on the emergence 
of a broad stationary out-of-equilibrium velocity distribution. Concluding 
remarks are drawn in Sec. [61 

2. The thermoelectric figure of merit ZT 

Let us focus our attention on a conductor in which both electric and heat 
current flow in one dimension (say, parallel to the x-direction). Assuming 
local equilibrium, a local entropy (per unit volume) s can be defined, and 
the rate of entropy production reads [lO| 

^ = J«5.(^)+Jp5.(-|), (1) 

in which J^ and Jp are the energy and particle current densities (fluxes) and 
9a;(l/T), —dx{iJi/T) the associated generalized forces (affinities), where T is 
the temperature and yU. the electrochemical potential. 

Assuming that the generalized forces are small, the relationship between 
fluxes and forces is linear and described b y th e phenomenological non equi- 



librium thermodynamic kinetic equations |10| . Ill 



Ju = L.uudAl^+Upd.A~^), (2) 



rji I ' fJ'l-' -^ \ rri 

Jp = ^p«5,.(^i^+L,,9,.(-|), (3) 

with La,i3 {a, (3 G {u, p}) Onsager coefficients. In the absence of magnetic 
fields, due to microscopic reversibility of the dynamics, the Onsager reci- 
procity relation L„p = Lpu holds. 

In analogy with the relation dQ = TdS, the heat current density Jq can 
be defined by the relation 

Jg = TJs, (4) 



with 

current density of entropy, and therefore 

Jq = Ju- jJ'Jp- (6) 

The entropy production rate equation can then be written in terms of the 
fluxes Jq and Jp and of the corresponding generahzed forces dx{l/T) and 
-{l/T)d,fi: 

s = JA[f)+J.{-^), (7) 

while the linear relationship between fluxes and forces reads as follows: 

J, = L,A (^) + L,p (-^) , (8) 

J. = LpA(f)+Lpp(-^), (9) 
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with Lpp = Lpp, Lgp = Lpg (Onsager relation), Lgp = L^p - fJ^Lpp, L 
Luu — '^fi'Lup + fi'^Lpp. Note that, if we call L and L the 2x2 Onsager 
matrices with matrix elements Lap {a,f3 G {u, p}) and L^s {l,^ ^ {lyP})y it 
turns out that detL = detL. 

The Onsager coefficients can be expressed in terms of more familiar quan- 
tities, the electric conductivity a, the thermal conductivity k, and the See- 
beck coefficient (thermopower) S. Let us first consider the case in which the 
thermal gradient vanishes, dxT = 0, and the system is homogeneous, so that 
the chemical potential pc is uniform. Since the electrochemical potential // is 
composed of a chemical part fi^ and an electric part /ig, /i = /ic + /ie, it turns 
out that for a homogeneous isothermal system dx/J, = d^He- The electric 
current Je = eJp, with e charge of the conducting particles, is then given by 
Je = cr£ = —{a/e)dxfie, with £ external electric field applied to the system. 
The quantities pc and /Xe cannot be determined separately by the theory of 



irreversible thermodynamics [12|: only their combination p = jjc + fJ^e ap- 
pears in the kinetic equations (|2]) and (|3]). Based on this equivalence, we can 
write Jp = Je/e = —{a/e^)dxfi even when /Xc 7^ 0, provided dxT = 0, whence 
Eq. (E]) gives 

a = —Lpp = —Lpp. (10) 



The heat conductivity n is defined as the heat current density per unit 
temperature gradient for zero electric current: Jq = —ndxT, at Jg = 0. 
Solving the two kinetic equations (|8]) and iQ simultaneously, we obtain 

1 det L 1 det L , , 

Finally, the Seebeck coefficient S is defined as the change in electro- 
chemical potential per unit charge, — 9^.yu/e, per unit change in temperature 
difference: S = —{l/e)dxii/dxT, at Jg = 0. We then obtain from Eq. ^ 

It is of course possible to eliminate the three Onsager coefficients Lqq, Lqp, 
and Lpp from the kinetic equations ([H]) and (jH]), and rewrite such equations 
is terms of the conductivities a and k, and of the thermopower S: 

Jq = -{K + TaS^)d,T - ^d^fi, (13) 

Jp = -^9x/i-— a.T. (14) 

e"^ e 

By eliminating d^fi from the above two equations one can express Jq in 
terms of Jp and d^T. It is then easy to derive an interesting expression for 
the entropy current density J^ = Jq/T J10l|: 

Js = eSJp - ^d,T, (15) 

from which the Seebeck coefficient can be understood as the entropy trans- 
ported (per unit charge) by the electron fiow. The second contribution to 
the entropy fiow, namely the term —{K/T)dxT, is independent of the particle 
current. 

The thermoelectric efficiency rj, of converting the input heat into output 
work, is determined by the non-dimensional figure of merit 

ZT = -— T. (16) 

To derive the relation between t] and ZT, we consider a one-dimensional 
system whose left/right ends are connected with left/right thermochemical 
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reservoirs, with small temperature difference AT = Tr — T^ and electro- 
chemical potential difference A/x = fiR — fiL- The efficiency rj is given, under 
steady state conditions, by the ratio of the time derivatives of the extracted 
work over the heat leaving the hot reservoir: 

Using Eqs. ( 1T^ and ( 1T^ to eliminate d^fi and Jg, we obtain 

where r^c- = 1 — Tr/Tl is the Carnot efficiency (here we assume Tl > Tr). 
The maximum efficiency for a given AT is derived after optimizing (fT8|) with 
respect to Jg: 

"-" = "° v/zr+i + r *"' 

The Carnot efficiency is therefore achieved in the limt ZT — )■ oo. 

Using Eqs. flTOj) . flTTl) . and flT2l) . we can express ZT in terms of the Onsager 
coefficients: 

ZT = = ^ P ~" ^^ pp' ('20') 

detL detL ■ 

The only thermodynamic restrictions to the Onsager coefficients come from 

the positivity of the entropy production, s > 0, which is a quadratic form in 

the generalized forces d.x{l/T) and —dx{fi/T) (see Eqs. ([I])-([3])) or 9^.(1/T) 

and —{1/T)dxfi (see Eqs. (I7])-([9])). Condition s > implies L^u^Lpp > 0, 

detL > in the first case, Lqq,Lpp > 0, detL > in the latter. Thus, the 

only restriction to the thermoelectric figure of merit is ZT > 0, so that in 

principle Carnot efficiency can be achieved. 

It is clear from Eq. ( l20l) that ZT diverges iff the Onsager matrix L 

(or, equivalently, L) is ill-conditioned, that is, when the condition number 

Ai(L)/A2(L) diverges, where Ai(L) and A2(L) are the largest and the smallest 

eigenvalue of L, respectively. The condition number diverges iff the quantity 

diverges. In this case the system ([2])- ([3]) (or, equivalently, the system ([8])- ([9])) 
becomes singular, and therefore J„ oc Jp. In short, the Carnot efficiency is 
obtained iff the energy and particle currents are proportional. 



3. Modeling thermochemical baths 

We consider a one-dimensional system whose ends are in contact with 
left/right baths (reservoirs), which are able to exchange energy and particles 
with the system, at fixed temperature T^ and electrochemical potential fia, 
where a = L,R denotes the left/right bath. 

The thermochemical reservoirs are modeled as infinite one-dimensional 
ideal gases. Therefore, particle velocities in the reservoirs are described by 
the Maxwell-Boltzmann distribution, 



where ks is the Boltzmann constant and m the mass of the particles. We use 
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a stochastic model of the thermochemical baths [13|: Whenever a particle 
of the system crosses the boundary which separates the system from the left 
or right reservoir, it is removed. On the other hand, particles are injected 
into the system from the boundaries, with rates 7q,. The injection rate 7q is 
computed by counting how many particle from reservoir a cross the reservoir- 
system boundary per unit time. That is to say. 



la = Pa dvvfo,{v)= Po,J^^-^, (23) 

Jo V 27rm 

with pa density of the ideal gas in reservoir a. Therefore, particles are in- 
jected into the system with velocity distribution 

where 9a{v) are step functions: 9l{v) = 1 if t' > 0, otherwise; 9ii{v) = 1 
ii V < 0, otherwise. We assume that injections from a macroscopic reser- 
voir are independent events and that the time interval between subsequent 
injections satisfies the Poissonian distribution, 

P„(t)=7„exp(-7„t), (25) 

so that the average time between injections is 1/ 'Ja- 
in order to relate the density pa to the electrochemical potential //„, it is 
convenient to write the grand partition function 

N 

(26) 



1 f A 

^ — < — e^"''" / dvmexp 

N=0 ' ^ 



-/3a ( -mv^ 



with A and A^ size and number of particles of the reservoir, respectively |j, 
Pa = l/lksTa) and h the Planck's constant. We then compute the average 
number of particles as 

so that 

P« - ^^ ^ • (28) 

Therefore, we can express the electrochemical potentials of the bath in terms 
of the injection rates: 

fj,a = kBTaln{\T^Pa), (29) 

with 

/?, 

\/2'KmkBTa 

de Broglie thermal wave length. Note that this relation, even though derived 

from the grand partition function of a classical ideal gas, can only be justified 

if particles are considered as indistinguishable. The 1/N\ term in the grand 

partition function (126|) is rooted in the above indistinguishability, of purely 



quantum mechanical origin [IJ]. The stochastic thermochemical baths used 
in our numerical simulations are based on Eqs. (!23|) . ( 124|) . (1251) . and (!29|) . 
The electrochemical potential /Xq, and the temperature T^ can be controlled 
by varying the injection rate 7^ and the temperature T^. 

4. One-dimensional non-interacting classical gas 

Let us first consider the simplest case of a one-dimensional gas of non- 
interacting particles. Assuming that also the reservoirs are one-dimensional 
and that the left/right contacts between system and reservoirs are identical 
and described as in Sec. [3l the particle current Jp reads 

Jp = lL I deuL{e)T{e) - 'jR deun{e)T{e), (31) 

Jo Jo 

where Ua(e) is the energy distribution of the particles injected from reservoir 
a and T(e) is the transmission probability for a particle with energy e to 



^It is of course understood that A is macroscopically large and that the thermodynamic 
hmit is eventuahy taken for the reservoir 
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transit from one end to the other end of the system, < T(e) < 1. Using 
Eq. ([M]), we obtain 

Uo.{e) = l3^e-^-\ (32) 

Furthermore, from Eqs. f l23|) and f l28|) we have 

7. = ^e"-. (33) 

After substitution of ( ^2]) and fl5^ into flMl) . we arrive to the following ex- 
pression for the particle current: 

1 /"°° 
Jp = j- de (e-^^(^-^^) - e-'^^^^-^^)) r(e). (34) 

Similarly, we obtain the heat currents Jq^a = Ju — fJ'aJp at the left and right 
reservoirs: 

J<i,c. = -^j deie-fi^){e-^-^'-^-^-e-M^-''-^)rie). (35) 

The thermoelectric efficiency is then given by (we assume Tl > Tr, fJ,R > fJ,L 
and consider only T(e) functions such that Jp> and Jgx > 0) 

^ J,,L - J,,R ^ im - /JL) ir de {e-M^-^-^ _ e-/^^(^-M.)) r(6) 
"^ J,,L /o"rfe(e-/ii)(e-'5i(-A^^)-e-/5«(-^«))r(e) ' ^ ^ 

When the transmission is possible only within a tiny energy window around 
e = e^,, the efficiency reads 

V = ^^^^^. (37) 

In the limit Jp — )■ 0, corresponding to reversible transport [6j, we get e* from 
Eq. dMD: 

'* = &-& ■ ''*' 

Substituting such e* in Eq. (l37|) . we obtain the Carnot efficiency rj = rjc = 
1 —Tr/Tl. Such delta-like energy-filtering mechanism for increasing thermo- 
electric efficiency has been pointed out in Refs. p, |6[. 

In the linear response regime, using a delta-like energy filtering, T(e) = 1 
in a tiny interval of width 5e around some energy e, otherwise, we obtain 

L6^(^6) .(,_^^ _ i:6(^6) g(,-_^^ _L(5e2^^(,_^) 

-"^un , , "^ ) -'^up -'^pM , , "^ ) pp ; ; ^ ' 

(39) 
10 
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Figure 1: Schematic picture of the model considered in our numerical simulations. 

where L is the length of system. From these relations we immediately de- 
rive that the Onsager matrix is ill-conditioned and therefore ZT = 00 and 
■J] = Tjc- We point out that the parameters e and 5e characterizing the trans- 
mission window, appear in the Onsager matrix elements ( 15^ and therefore 
are assumed to be independent of the applied temperature and electrochem- 
ical potential gradients. On the other hand, the energy e^ in Eqs. ( 137|) .( 138|1 
depends on the applied gradients. There is of course no contradiction since 
( 137|) ■ ( 138|) have general validity beyond the linear response regime. 

5. One-dimensional interacting classical gas 

Let us now turn to the interacting case. We consider a one-dimensional, 
di-atomic disordered chain, of hard-point elastic particles with coordinates 
Xi e [0,L], L being the system size, velocities Vi and masses rrii G {m, M } 
randomly distributed. The particles interact among themselves through elas- 
tic collisions only. A schematic picture of the model is drawn in Fig. [T] Since 
we are considering a purely mechanical model, strictly speaking we are going 
to investigate thermodiffusion rather than thermoelectricity. On the other 
hand, we assume that the particles are charged and that the Coulomb repul- 
sion is screened and modeled by a short-range hard-core interaction (elastic 
collisions). Therefore, our model is relevant also for thermoelectricity. Nu- 
merical results obtained in Ref. |9| suggest that, for mass ratio M/m 7^ 1, 
the figure of merit ZT diverges in the thermodynamic limit. □ 

Let £ be a reference unit length which we take 1 in simulations. In our 
numerical simulations we set // = (/i^ + /i_R)/2 = 0.2 [h? /m(?] and T = 
{Tl + Tr)/2 = 3.0 [h? /m(?kB\, and consider /i^,, /i^j, T^, Tr slightly different 
from the mean values /i, T to drive finite currents Jp and J„. We assume that 



^The two masses must be different in order to have ergodic and mixing dynamics, so 
that thermalization within the system occurs. For equal masses the dynamics is integrable 



andZT= 1 0. 
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the mass of each particle injected by the left or right bath is chosen randomly 
and with equal a priori probabilities between the two possible values m and 
M. The average currents J„ and Jp are computed at the contacts between 
system and baths: If, in a period of time t the left bath injects Ni particles 
with masses rrij and velocities Vj , j = l,...,iVj, and absorbs Na particles 

with masses rrij and velocities w , j = l,...,Na, then in the large t limit 
the currents Jp and J„ are given by 

J, = j(N, - iVJ, (40) 

\j=i j=i / 

Note that in the steady state, due to particle and energy conservation, these 
currents are equal to the corresponding currents computed for the right bath. 
Then the Onsager matrix elements from which a, S, k, and ZT can be readily 
derived, are obtained from Eqs. ([8]) and ([9]). We set the mass ration M/m = tt 
and calculate currents up to L = 80[i], corresponding to an average number 
of particles inside the system (A^) f« 515. 

In Fig. [2], we present our numerical results for the transport coefficients. 
A power law dependence for a/L, n/L, and ZT is observed |f|. In particular, 
the figure of merit ZT increases with increasing the systems size, ZT oc 
L°-^^. Correspondingly, the condition number [Tr(L)]^/det(L) (see Eq. (12T|) ) 
diverges, as expected from the general theoretical considerations of Sec. [2l 

These numerical results naturally raise a question: Is the mechanism 
leading to high ZT quality factor for interacting gases related to the delta- 
like mechanism p, |6| shortly discussed in Sec. H] for the non-interacting ideal 
gas? To address this question, we measure the particle current at the position 
X G [0, L] as 

/■oo 

Jp= / dED{E), (42) 

DiE)=DLiE)-DR{E), (43) 

where the "transmission function" Dl[E) is the density of particles with 
energy E crossing x and coming from the left side, while Dr{E) is the density 



^Numerical data are consistent with those reported in Ref. [9| for the same mass ratio. 
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Figure 2: Thermoelectric transport properties. The quantities <j/L, k/L, ZT, and ^^^^1 
show a power law dependence on the system size L. 
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Figure 3: D{E) calculated for L = 10,40 and 80. 

of particles with energy E from the right side. We inquire how D{E) changes 
as a function of L, in particular if D{E) becomes more and more delta-like 
(peaked in energy) when increasing L. The transmission function D{E) is 
shown in Fig. [3l at x = L/2 and for different system sizes Q There is no sign 
of narrowing of D{E) when increasing the system size. We can therefore 
conclude that the mechanism leading to the large ZT values observed in 
Fig. [2] must be different from the energy filtering discussed in Refs. p, l6| . 

To understand the mechanism, we first write the particle and energy 
currents as 



Jp = v{x,t)p{x,t), Ju 



-mv{x,t)^p{x,t), 



(44) 



where x G [0,L], the overbar denotes time-averaging, and v{x,t), p{x,t) are 
respectively the particles velocity and density at the position x and time t. 
If the relaxation time scales for density and velocity are well separated, then 
expressions (jUj) can be approximated as: 



Jp ^ v{x,t) X p{x,t), Ju ^ -mv{x,t)^ X p[x,t). (45) 

In our model, this is satisfied. For instance, in the case of {pl.pr) = 



^Note that, while Jp is position-independent due to conservation of particles, D{E) 
depends on x. However, we have checked that similar behaviors of D{E) are obtained for 
different values of x. 
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Figure 4: F{v) calculated for the parameters {fiL,iJ,B.) = (0.24,0.16) [/i^/mi!^],T = 
3.0 [/i^/mf^fcs] with the system size L = 40. The solid line corresponds to the the numer- 
ical data, at a; = L/2. The mean velocity is v(L/2,t) ^ 0.010 [h/m£]. The dashed line is 
the modified Maxwcll-Boltzmann distribution FMiv) which fits F{v) with the parameters 
(v) = v{L/2,t) = 0.010[h/mf] and i' = 1.15[h/mi]. For these parameters, Eq.gZl) yields 
(v^) ^ 0.041, which is comparable to the exact numerical data (v^) ~ 0.047. The inset 
shows the behavior of F{v + (v)) and F{—v + (v)) in the semi-log scale. These two curves 
completely overlap each other, as expected from Eq. pS)) . Note that the tails deviate from 
the Gaussian form. 

(0.24,0.16) [h'^/mi%T = 3. [/iVrnf ^^/gg], and L = 40 [£], we get at x = L/2 
the time-averaged velocity v{L/2,t) ~ 0.010 [h/mi], and the time-averaged 
density p(L/2,t) ^ 6.43 [l/£], while v{L/2,t) p{L/2,t) ^ 0.0641 [h/mi^]. 

From the discussion of Sec. [2l it is clear that ZT diverges when J„ oc Jp. 
According to Eq. fHSjl . this is the case when v^ oc v. Since we are interested 
in the steady-state transport properties and we are considering systems with 
the mixing property, it is natural to assume that the time-averages f " equal 
the ensemble averages (w") = j_ dvv"'F{v), with F{v) velocity distribution 
function for the steady state. At equilibrium (T^ = Tr, /i^ = /ir), the system 
thermalizes and F{v) is the Maxwell-Boltzmann distribution (122|) at any x. 
In the linear response regime, we assume that F{v) is given by a "modified 
Maxwell-Boltzmann distribution" , 



Fm{v) 



m 



27ikBT 



exp 



'm*{v — {v)y 
2kBT 



(46) 



where the mean velocity (v) and the effective mass m* are fitting parameters, 
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and T ^ Tl ^ Tfj. That is to say, we assume that the out-of-equihbrium 
stationary distribution (1461) differs from the equihbrium Maxwell-Boltzmann 
distribution only in the position (v) of the peak, while the Gaussian shape is 
unchanged. As shown in Fig. HI such assumption is in good agreement with 
the numerically computed F{v) close to the peak of the distribution, while 
the tails show deviations from fH^ . Nevertheless, such deviations do not 
affect too much the values of (v) and (v^) and Eq. (1461) is very convenient for 
analytical considerations and to unveil the mechanism at the origin of the 
large thermoelectric efficiencies observed in the hard-point gas model. 
From Eq. (146|) we obtain 

{v^) = {vy' + 3u^{v), (47) 

where 

-/? («) 

is the width of distribution (l46l) . We obtain (v^) oc (v) when u ^ (v), 
that is, in the broad-distribution limit. It is clear from Fig. Hlthat, for the 
one-dimensional interacting hard-point gas, indeed u ^ {v)-u 

6. Conclusions 

We have studied numerically the coupled particle and energy transport in 
a prototype model of interacting one-dimensional gas: the disordered, hard- 
point gas. There is numerical evidence that the ZT quality factor diverges 
with increasing the system size. We explain this result in terms of the emer- 
gence of a broad velocity distribution of the particles transmitted across the 
sample. This mechanism first of all requires local equilibrium, which is ex- 
pected to take place in systems with the mixing property. We also make a 
couple of assumptions which are quite natural in many-body systems: the 
separation of the relaxation time scales of density and velocity in Eq. fj^5|) . 
and the modified Maxwell-Boltzmann form of the velocity distribution P6|l . 
On the other hand, since ZT = [aS"^ / k)T and Fig. [2] shows that the Seebeck 



^ Note that the delta- hke limit of Eq. (pS)) . i^ <C (w), is incompatible with the linear 
response regime plus approximation (j45p since, if Jp ex {v) is a linear function of the 
applied temperature and electrochemical potential gradients, the same cannot hold for 
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coefficient is practically constant, the anomalous behavior of a and k [15 
is crucial to obtain a diverging ZT. The relationship between the broad 
velocity-distribution mechanism and the anomalous behavior of the trans- 
port coefficients must be clarified. In particular, further investigations are 
required to understand whether this mechanism could be applied to systems 
with the mixing property but without anomalous transport. It might indeed 
be possible to find systems in which a, n and ZT eventually converge to finite 
but large values, when increasing the system size. Therefore, our mechanism 
could be also relevant in more realistic interacting systems with the mixing 
property. 
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